Measurement of Scattering Rate and Minimum Conductivity in Graphene 
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The conductivity of graphene samples with various levels of disorder is investigated for a set of 
specimens with mobility in the range of 1 — 20 x 10 3 cm 2 /V sec. Comparing the experimental 
data with the theoretical transport calculations based on charged impurity scattering, we estimate 
that the impurity concentration in the samples varies from 2 — 15 x 10 11 cm -2 . In the low carrier 
density limit, the conductivity exhibits values in the range of 2 - 12e 2 /h, which can be related to 
the residual density induced by the inhomogeneous charge distribution in the samples. The shape 
of the conductivity curves indicates that high mobility samples contain some short range disorder 
whereas low mobility samples are dominated by long range scatterers. 
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The recent discovery of graphene [l|, a single atomic 
sheet of graphite, has created an intense research effort 
on this new material. The uniqueness of the electronic 
band structure of graphene is exemplified by the charge 
neutral Dirac points located at the Brillouin zone cor- 
ner. The linear dispersion relation of the bands cross- 
ing these Dirac point can be interpreted in analogy to 
relativistic fermions of vanishing mass. The nontrivial 
Berry phase associated with the graphene band struc- 
ture accounts for unusual half integer shifted quantum 
Hall effect @, H in the presence of a magnetic field. 
This chiral nature of the Dirac fermions is also expected 
to induce nontrivial charge transport phenomena, such 

as a minimum conductivity @, 11 S, 0, 0] 

and weak 



anti- localization 0, 0, [l(J 11 1. Indeed, transport mea- 
surements indicate the existence of such a minimum con- 
ductivity [2, 4] and the suppression of localization [13] in 
some graphene samples. Despite qualitative agreement, 
however, these experimental observations do not agree 
quantitatively with the theoretical expectations @, U Ell- 
in order to understand these and other low energy trans- 
port phenomena in graphene, it is essential to study rel- 
evant carrier scattering mechanisms. 

In this Letter, we report low magnetic field transport 
measurements on 19 graphene devices with various levels 
of disorder. We infer the scattering mechanisms using 
the density dependence of the mean free paths and phase 
coherence lengths and comparing with the corresponding 
transport theory. The sample mean free path is extracted 
from the conductivity measurements, and we find that at 
high carrier density, two different scattering mechanisms 
determine the density-dependent conductivity. 

The graphene pieces used in this work are extracted 
from bulk graphite crystals and deposited onto SiC^/Si 
substrate u sing the mechanical method described in pre- 
vious work [13[ . Typically, graphene pieces of lateral size 
3-20 /im are chosen for device fabrication. For a few 
large flake samples (lateral size > 20 /im), we explic- 
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FIG. 1: (color online) (a) Resistivity of four representative 
graphene samples as a function of applied gate voltage. All 
data were measured at temperature 1.6 K. (b-d) Optical mi- 
croscope images of typical devices. The scale bar represents, 
1, and 10 ^m, for (b) and (c, d), respectively 



itly fabricated well-defined two to three Hall-bar geome- 
try devices having the same width but varying lengths. 
These multiple devices from the same graphene piece pro- 
vide means to check the consistency of our conductiv- 
ity measurements. The metallic electrodes are defined 
on the sample using electron beam lithography followed 
by Au/Cr (30/3 nm) evaporation and a lift-off process. 
The degenerately doped silicon substrate serves as a gate 
electrode with 300 nm thermally grown silicon oxide act- 
ing as the gate dielectric. By applying a gate bias volt- 
age, V g , the charge density of the sample can be tuned. 
The conductivity of a sample is measured either in the 
Hall bar geometry considering the width and length of 
the samples in between the voltage probes (Fig. 1(b)) 
or by van der Pauw method for samples with irregular 
shape (Fig. 1(c)). We also use Hall-bar shaped devices 
(Fig. 1(d)) in order to minimize systematic measurement 
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FIG. 2: (color online) Mobility estimated by applying Drude 
model to the data in Fig. 1. The graphs shown by broken 
lines are the region where simple Drude model fails to hold. 
The inset shows the scattering time in this model. The solid 
lines are for electrons and dashed lines are for holes. 



errors. All samples are measured using lock-in amplifiers 
at an excitation current less than 50 nA to minimize any 
heating effect on the samples, which are kept in the he- 
lium vapor of a cryostat. 

Fig. 1(a) shows the resistivities, p, of five representa- 
tive samples as a function of V g . The resistivity curves 
are largely symmetric around a particular gate voltage 
V g = Vuirac and show a maximum at this value. The 
finite value of Vuirac indicates that there exists an un- 
intentional doping of the graphene samples [14j | whose 
origin may be electrostatic or caused by charged impu- 
rities. The actual carrier density n in graphene induced 
by the gate voltage in the presence of impurity doping 



is then obtained from n = C g (V g — Voiracj/s. where the 
gate capacitance C g = 115 aF//im 2 is deduced from a 
separate Hall measurement. We observe that p decreases 
monotonically as \n\ increases, however the sharpness of 
the dips at n — varies drastically from sample to sam- 
ple. 

In order to analyze the difference between samples 
quantitatively, we employ the semi-classical Drude model 
to estimate the mobility of samples, p = (enp)" 1 . Fig. 2 
display p as a function of n obtained from Fig. 1(a). Al- 
though p defined in this way diverges at V g = Voirac (the 
portion of curves represented by broken lines) due to the 
trivial reason that n = nominally at the Dirac point, 
the limiting value of the mobility, p^, defined in the large 



density limit (n ~ 4 x 10 12 cm -2 ) serves as a useful way 
to characterize the sample quality. We have measured a 
total of 19 samples in this experiment where ph ranges 
2, 000 — 20, 000 cm 2 /Vsec. Quite generally, for samples of 
poorer quality, the Dirac points shift more from V g = 0, 
indicating a larger, unintentional charge doping of un- 
known origin [15j. The inset to Fig. 2 shows explicitly 
the scattering time r of the samples shown in Fig. 1, with 
r = pEp/(ev F ) by definition. The corresponding mean 
free path of the samples can be obtained from I = vpr, 
where the Fermi velocity vf ~ 10 6 m/sec [4|. We note 
that the resulting mean free path for both electrons and 
holes strongly depends on the density n, ranging from 10- 
500 nm in the most of the samples and density ranges. 

The quality of the samples can be further differenti- 
ated by comparing their behavior of carrier density de- 
pendent conductivity. Fig. 3 shows a — p" 1 as a func- 
tion of V g . Evidently, the samples with poorer quality 
(pL < 5,000 cm 2 /Vsec) show a very broad and smooth 
maximum in p(V g ) near the Dirac point followed by a 
linear relationship of ^(V^) in the large density limit. On 
the other hand, when the mobility of a sample is larger 
than ~ 10,000 cm 2 /Vs, <r{V g ) forms a cusp around the 
Dirac point followed by a sublinear increase both in the 
electron and hole regimes. Since the scattering time r 
goes as t ~ cm -1 / 2 , this difference in the density depen- 
dence of a in the two different mobility regimes strongly 
indicates that different scattering mechanisms might be 
dominating the high carrier density transport in samples 
from the two limiting groups. Recent numerical calcula- 



tions based on Ref. [7], [14| indeed show that uin) changes 



from a linear dependence to a sublinear dependence as 
the scattering mechanism changes to short range scatter- 
ing, such as atomic defects in the lattice, from long range 
scattering, such as ionized impurity scattering. 

Comparing our experimental data in a quantitative 
fashion with these theoretical calculations TJ], we es- 
timate that the charged impurity concentration varies in 
the range of 2 — 15 x 10 11 cm" 2 . Compared to recent Ra- 
man data on graphene samples, where the typical impu- 
rity density rii is estimated to be n; ~ 3 x 10 11 cm" 2 , this 
value appears to be quite reasonable [llij . We note that 
the relevant scattering times entering Raman scattering 
level broadening and dc conductivity are not necessarily 
identical, and typically the level broadening is larger than 
the transport broadening for the same impurity density. 
This explains why the Raman impurity density estimate 
falls toward the lower end of our transport estimate. 

Simple theoretical arguments show that long ranged 
charged impurity scattering is much more effective for 
scattering of carriers than short-range scattering 0, [3] • 
For example, given the same impurity density rii ~ 
5 x 10 11 cm" 2 , I > 1 pm for short-range scatterers but 
I ~ 50 nm for charged impurity scatterers. Since the size 
of our samples L is on the order of several microns, trans- 
port can be considered to be in the diffusive regime. A 
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FIG. 3: (color online)The conductivity of five representative 
samples shown in Fig. I as a function of applied gate volt- 
age (solid lines). For clarity, the each curves are displaced. 
The horizontal dashed lines indicate the zero conductance for 
each curves. Dotted curves are the corresponding theoretical 
results as explained in the text. The only fit parameter in the 
theory is the density of the random charged impurity centers 
in the substrate, which is taken to be (top to bottom in units 



of 10 1 



2.2, 4.0, 4.6, 9.7, 14.5. The inset shows the 



detailed view of the density dependent conductivity near the 
Dirac point for the data in the main panel. 



Boltzmann transport theory with charged impurity scat- 
tering can then be applied to solve for the conductivity 
self-consistently, and we obtain [l7| 



h rii 



C^-— for n < 71* 

n rii 



(i) 



where the carrier density is related to the gate voltage 
by n — CgVg/e + n; n* describes the self-consistent resid- 
ual carrier density considering electron and hole pud- 
dle formation induced by the charged impurities; fi is 
the induced carrier density in graphene by the charged 
impurity; and C is a dimensionless numerical parame- 
ter describing the strength of scattering by the potential 
fluctuation considering full random phase approximation 
(RPA) screening. Considering the dielectric screening 
from the SiC>2 substrate underneath the graphene sam- 
ples, we estimate C w 20 17]. Within the RPA, the 
values of n* and n can be self-consistently determined by 
the charged impurity concentration rii, and d, the dis- 
tance between the graphene sheet and the charged impu- 
rity layer [T3, - 



We emphasize that in this model with a single free pa- 
rameter rii, Eq. [T] describes: (i) the linear tail mobility, 
(ii) the offset voltage Voirac = —en/C g , (hi) the mini- 
mum conductivity given by m 20e 2 /h(n*/n l ) and (iv) the 
plateau width i.e. the range in gate voltage over which 
the conductivity saturates at its minimum value. Since 
these features are determined by the impurity concentra- 
tion only, we now understand qualitatively that the be- 
havior of conductivity of the cleaner samples (i.e. those 
with lower charged impurity density), in addition to hav- 
ing a higher slope, also have narrower minimum conduc- 
tivity plateaus, smaller values of Voirac, and higher val- 
ues of the minimum conductivity (i.e., o{V g — Voirac))- 
On the other hand, the dirtier samples have lower mobili- 
ties, wider plateau widths, and larger gate voltage offsets. 

For a quantitative comparison, the theoretical model 
(Eq. [1]) is now compared to the experimental data 
(Fig. 3). To avoid having too many fitting parameters 
in the theory, we just set d w 1 nm (which is typical for 
SiC>2 substrate [l7(), keeping only rii as a single fitting 
parameter. This yields corresponding charged impurity 
concentrations for the samples in Fig. 3 ranging between 
2 — 15 x 10 11 cm" 2 , where a reasonable agreement be- 
tween experiment (solid curves) and theoretical curves 
can be seen. We note that that the some discrepancies 
between the experimental data and the simple model can 
be further understood: for example the sub-linear be- 
havior seen in high quality samples (K17 and K12) is 
consistent with short-range scatte ring that becomes im- 
portant only at very high density 14|. The electron and 



hole asymmetry often observed in poor quality samples 
(K130) can be explained by drift of charged impurity in 
the substrate. The large gate voltages used in the ex- 
periment (V g ~ 50 V) can shift the impurity trapped in 
the substrate, resulting in the adjustment of d for elec- 
trons and holes. Had we introduced two additional fit pa- 
rameters (one to parameterize the amount of short-range 
scattering and another one to parameterize d the distance 
between the charged impurities and the graphene sheet) 
we would have obtained exact quantitative agreement. 
Given the reasonable agreement between theory and ex- 
periments for several different samples in Fig. 3, we do 
not see any advantages in obtaining an exact quantitative 
agreement by doing this additional data fitting. 

We now address the minimum value of the conductiv- 
ity o"mm observed in our samples. Fig. 4 displays a m i n 
versus the high-density, limiting mobility value (i^ for 
all our samples. The measured values of a m in he in the 
range of 2 — 12e 2 /h. Most samples have a m in > 4e 2 //i, 
in contrast to an earlier study [2| reporting a clustering 
of cr m i n ~ 4e 2 /h (horizontal dotted line). Several theo- 
ries 0, 0] based on homogeneous graphene sheets predict 
a universal value of 4e 2 /irh for a m i n (horizontal dashed 
line), which is almost an order of magnitude lower than 
in most of our samples. The sample dependent a m i n was 
recently argued to arise from the inhomogeneous charge 
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FIG. 4: The minimum conductivity and limiting mobility val- 
ues at metallic regime for all 19 samples measured at 1.6 K. 
Filled symbols indicate two devices, K145 and K146 made out 
of the same graphene flake as shown in the optical microscope 
image in the inset. 



distribution near the Dirac point [i, 14 1. More recently, 
theoretical work based on numerical calculation also in- 
dicates a non-invariant cr m ,- n scaling with the system size 
in clean limit [l|| [2(J. Further careful consideration 
is needed for disordered samples. In the low density 
limit near the Dirac point where the carrier concentra- 
tion becomes smaller than the charged impurity density, 
the system breaks up into puddles of electrons and holes 
where a duality in two dimensions guarantees that locally 
transport occurs either through the hole channel or the 
electron channel. Percolation resistance of these electron 
and hole puddles are the key to understanding the ex- 



perimentally observed minimum conductivity 21]. The 
charged impurities induce a density distribution in the 
graphene sample where the presence of both electrons 
and holes carriers implies that both positive and nega- 
tive fluctuations are screened by graphene 17] . It is the 
residual density given by n* that is responsible for the 
minimum conductivity, and therefore we expect that the 
value of the minimum conductivity will depend on the 
impurity concentration. Very recent scanning single elec- 
tron transistor microscopy demonstrated such electron- 
hole puddle formation near the charge neutral point in 
graphene [221 ]. 

In conclusion, we report zero-field transport properties 
of more than a dozen graphene samples with different 
mobility values. The minimum conductivity is found to 
be strongly sample dependent, yielding values an order 
of e 2 /h with a non-universal prefactor. Density inhomo- 
geneities across the graphene sheet may be responsible 
for this large spread in conductivity. 
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